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EXTREMAL COVARIANT POVM'S 



GIULIO CHIRIBELLA AND GIACOMO MAURO D'ARIANO 
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^Nj , Abstract. We consider the convex set of positive operator valued measures 

(POVM) which are covariant under a finite dimensional unitary projective 
representation of a group. We derive a general characterization for the ex- 
tremal points, and provide bounds for the ranks of the corresponding POVM 
_ ^ densities, also relating extremality to uniqueness and stability of optimized 

>««^ 1 measurements. Examples of applications are given. 

m 



1. INTRODUCTION 



An essential step in the design of the new quantum information technology [hi 
is to asses the ultimate precision limits achievable by quantum measurements in 
^SJ ■ extracting information from physical systems. For example, the security analysis 

^O ! of a quantum cryptographic protocol Q is based on the evaluation of the limits 

^P ' posed in principle by the quantum laws to any possible eavesdropping strategy. 

A general method to establish such limits is to optimize a quantum measurement 
according to a suitable criterion, and this is the general objective of the so-called 
quantum estimation theoryf^ ^ . Different criteria can be adopted for optimizing the 
measurement, the choice of a particular one depending on the particular problem 
at hand. Moreover, many different optimization problems often share the same 
jrt ' form, e. g. they resort to the maximization of a concave function on the set 

^ I of the possible measurements. We remind that measurements form a convex set, 

"sJ ' the convex combination corresponding to the random choice between two different 

apparatuses. Since a concave function attains its maximum in an extremal point, 
it is clear that the optimization problem is strictly connected to the problem of 
characterizing the extremal points of the convex set. 
j^ ' The quantum measurements interesting in most applications are covorianflij 

with respect to a group of physical transformations. In a purely statistical descrip- 
tion of a quantum measurement in terms of the outcome probability only — i. e. 
without considering the state-reduction — the measurement is completely described 
by a positive operator valued measure (POVM) on its probability space. In terms 
of POVM's, " group-covariant" means that there is an action of the transformation 
group on the probability space which maps events into events, in such a way that 
when the measured system is transformed according to a group transformation, 
the probability of a given event becomes the probability of the transformed event. 
Such scenario naturally occurs in the estimation of an unknown group transforma- 
tion performed on a known input state, e. g. in the estimation of the unknown 
unitary transformation pi o, in the measurement of a phase-shift in the radiation 
field 0, ItI , or in the estimation of rotations on a system of spins [B . A first tech- 
nique for characterizing extremal covariant POVM's and quantum operations has 
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been presented in Ref. [pi inspired by the method for characterizing extremal cor- 
relation matrices of Ref. [|0|, in particular, classification of extremal POVM's has 
been presented for the case of trivial stability group, i. e. when the only trans- 
formation which leaves the input state unchanged is the identity. Here we solve 
the characterization problem for extremal covariant POVM's in the general case of 
nontrivial stability group, providing a simple criterion for extremality in Theorem^ 
in terms of minimality of the support of the seed of the POVM, presenting iff condi- 
tions for extremality in Theorem 0, and providing bounds for the rank of extremal 
POVM's (in the following we will define the rank of a POVM as the rank of its 
respective density: sec Eq. (m for its definition). We show that, contrarily to the 
usual credo, the optimal covariant POVM can have rank larger than one. Indeed, 
there are group representations for which covariant POVM cannot have unit rank, 
since this would violate a general bound for the rank of the POVM in relation to 
dimensions and multiplicity of the invariant subspaces of the group. In the present 
paper we adopt the maximum likelihood optimality criterion, which, however, as 
we will show, is formally equivalent to the solution of the optimization problem in 
a very large class of optimality criteria. Other issues of practical interest that we 
address are the uniqueness and the stability of the optimal covariant POVM. The 
whole derivation is given for finite dimensional Hilbert spaces: as we will show in a 
simple example, it can be generalized to infinite dimensions, however, at the price 
of making the theory much more technical. 

The paper is organized as follows. After introducing covariant POVM's and their 
convex structure in Section |[ the main group theoretical tools that will be used for 
the characterization of covariant POVM's are presented in Section y. In Section |j 
we give a characterization of extremal covariant POVM's in finite dimension with 
a general stability group, deriving an algebraic extremality criterion, along with a 
general bound for the rank of the extremal POVM's in terms of the dimensions of 
the invariant subspaces of the group and of the stability subgroup. Properties of 
extremal POVM's in relation with optimization problems are analyzed in Section || 
where also the issues of uniqueness and stability of the optimal covariant POVM's 
are addressed. Finally, examples of application of the theory to estimation of ro- 
tation, state, phase-shift, etc. are given in Section @, providing extremal POVM's 
with a non trivial stability group and giving examples of optimization problems 
with solution consisting of extremal POVM with rank greater than one. 

2. Convex structure of covariant POVM's 

The general description of the statistics of a measurement is given in terms of a 
probability space X — the set of all possible measurement outcomes — equipped with 
a cr— algebra a{X) of subsets B C X and with a probability measure p on a(X). 
Each subset B G cr(X) describes the event "the outcome x belongs to B" and the 
statistics of the measurement is fully specified by the probability measure p, which 
associates to any event B its probability p(B). 

In quantum mechanics the probability p(B) is given by the Born rule 

(1) piB) - Tr[pP(B)] 

where p is a density operator (i.e. a positive semidefinite operator with unit trace) 
on the Hilbert space J^ of the measured system, representing its state, whereas P 
is the POVM of the apparatus, giving the probability measure p for every given 
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state p of the quantum system. Mathematically a POVM P : cr(X) -^ B{J^f) is a 
positive operator valued measure on cr(X), namely it satisfies the following defining 
properties 

(2) < P(B) <I VB e a{X) 

OO 

(3) P(U,^iB,)=^P(B,;) V{BJ disjoint 

(4) P(X) = /. 

Notice that the set of POVM's for a(X) is a convex set, namely, if Pi and P2 are 
POVM's for fT(X), then also APi + (1 - A)P2 is a POVM for ct(X) for any < A < 1. 
The measurement described by the POVM APi + (l — A)P2 corresponds to randomly 
choosing between two different measuring apparatuses described by the POVM's 
Pi and P2 respectively. The extremal points of such convex set of POVM's — the 
socalled extremal POVM's — correspond to measurements that cannot result from 
a random choice between different measuring apparatuses. 

In the following we will focus attention to the case of probability space X given 
by the quotient G/Gq of a compact Lie group G with respect to a subgroup Gq. 
Physically, this situation arises when the POVM is designed to estimate a state of 
the group-orbit {UgpU^ \ g G G} of a given state p, with the group G acting on 
the Hilbert space J^ of a quantum system via the unitary projective representation 
R(G) = {Ug I g G G}. In such case, in fact, the probability space of the POVM is 
exactly X = G/Gp, and Gq = {h e G \ UhpUl — p} is the stability group of p, 
whence the points of the orbit are in one to one correspondence with the elements 
of X = G/Gq. Notice that in the following the fact that the representation is 
projective is inconsequential, whence there will be no need of reminding it. 

An important class of measurements with X = G/Gq is described by the covari- 
ant POVM's M, namely those POVM's which enjoy the property 

(5) P(gB) = UgPiml VB e a(X), Vg G G, 

where gB> ^ {gx \ x E B}. Any POVM P in this class is absolutely continuous with 
respect to the measure da; induced on X by the normalized Haar measure dg on 
the group G, and admits an operator density M, namely 

(6) M : X^B{J^), P(B) = I AxM{x). 

Jb 

For a covariant POVM, the operator density has the form [Q 

(7) M{x) = Ug^,jEUl^^^, 

where g{x) G G is any element in the equivalence class x G X = G/Gq, and S is 
an Hermitian operator satisfying the constraints 

(8) S>0, / dgUgEUl =/ 

Jg 

(9) [S,C/,,]=0 V/igGo. 

The operator S is usually referred to as the seed of the covariant POVM[|ll[. 

Notice that the constraints (g|) are needed for positivity and normalization of the 
probability density, whereas identity (||) guarantees that M{x) = Ug{^x-^'EU ,. does 
not depend on the particular element g(x) in the equivalence class x. It is easy 
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to see that the constraints (g) and (^ still define a convex set C, namely, for any 
Si , S2 e C and for any < A < 1 one has ASi + (1 — A)S2 G C. Precisely, the convex 
set C is the intersection of the cone of positive semidefinite operators with the two 
affine hyperplanes given by identity (0) and by the normalization condition in Eq. 
(||). Since a covariant POVM is completely specified by its seed S as in Eq. (0), 
the classification of the the extremal covariant POVM's resorts to the classification 
of the extremal points in the convex set C. 

3. Group theoretic tools 

Let G be a compact Lie group, with invariant Haar measure dg normalized as 
L-, d (7 = 1, and consider a unitary representation R(G) = {Ug \ g G G} on a finite 
dimensional Hilbert space Jff. Then J^ is decomposed as direct sum of orthogonal 
irreducible subspaces as follows 

(10) ^ = 00^(A'), 

fies 1=1 

5 denoting the collection of equivalence classes of irreducible components of the 
representation, the classes being labeled by the Greek index fj,, whereas the Latin 
index i numbers equivalent representations in the same class. Let T^^ : J^- — )■ 
Jfff^ denote invariant isomorphisms connecting the irreducible representations of 
the equivalence class /i of dimension d^, namely for any i,j = 1, . . . ,771^ TJ^j : 
J^j -^ Jif^ is an invertible operator satisfying the identity 

(11) UgT^^^Ul = T^^, VgeG. 

Consistently with this notation T^^^ will denote the projection operator on ,J^ . 
Since all subspaces Jff^ are isomorphic, we can equivalently write 



(12) ^jqM^jf^(^^^ 



(m) _ 

■ ^ M ' 

1 



where c^ denotes the representation space, i.e. an abstract d^-dimensional sub- 
space where a representation of the class /i acts, while ^^ denotes the multiplicity 
space, i.e. a TTi^-dimensional space which is unaffected by the action of the group. 
In this way, the decomposition (lin) can be written in the Wedderburn's formjl3| 

(13) ^ = 0jr^®^^. 

Due to Schur lemmas, an operator O in the commutant of the representation 
R(G) can be decomposed as follows Q 

whereas, in terms of the decomposition (|3|) one has 
(15) O = ®^6s {I^. ® Op) ' 

/^ denoting the identity on the representation space J^, and O^ e S(^p) being 
a suitable set of operators on the multiplicity spaces ^p. 
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In this paper we will consider covariant POVM's with X = G/Go where both G 
and Go are compact Lie groups, represented on the Hilbert space M' by the unitary 
representations R(G) = {Ug | 5 G G} and R(Go) = {Uh \ h € Go}- We wiU denote 
with S and So the equivalence classes of irreducible representations of R(G) and 
R(Go) respectively. The constraints Mm) can be rewritten in a remarkably simple 
form using the decompositions of J^ in irreducible subspaces under the action of 
R(G) and R(Go). In fact, due to the invariance of the Haar measure dg, the integral 
in (ph belongs to the commutant of R(G). Rewriting the constraint (H) by using 
(p^, one get easily: 

(16) Tr[7;f)s]=d^%, VmgS, Vi, j = 1, . . . , m^. 

Moreover, according to (ph and (0) , the operator S must be a positive semidefinite 
operator in the commutant of R(Go) (|9|), then we have 

(17) E = (B.eSoiIu(^XlX,), 
where Xi, is an operator on the multiplicity subspace ^i,. 

4. Extremal covariant POVM's with a nontrivial stability group 

In this section we will classify the extremal points of the convex set C of covariant 
seeds, namely the convex set of operators that satisfy both conditions (||) and (||). 
For the characterization of the extremal points of a convex set we will use the well 
known method of perturbations. We will say that the operator is a "perturbation" 
of a given S g C if and only if there exists an e > such that S + t0 € C for any 
t S [— e, e]. With such definition one has that an operator S is extremal if and only 
if its unique perturbation is the trivial one, namely if is a perturbation of S then 
9 = 0. 

Let's start with a simple lemma which is useful for the characterization of the 
perturbations of a given seed S. 

Lemma 1. Let S G B{Jf) be a positive semidefinite operator. Then, for any 
Hermitian O G S(^) the condition 

(18) 3e>0: Vt G [-e, e] S + tO > 
is equivalent to 

(19) Supp(e) C Supp(S). 

Proof. Suppose that the condition ( pq ) holds. Then for any \(j)) G Ker(S) one 
necessarily has {(t)\<d\(t>) = 0. Therefore, for any vector |^) G ^ one has: 

mem = \m{E+tQm\ < -^y/mw+mWJm+teM = o. 

Hence Ker(S) C Ker(9), implying that Supp(6) C Supp(S). Conversely, suppose 
that (noh holds. Let's denote by A the smallest nonzero eigenvalue of S and by 1 10| | 
the norm of O, then condition ( p^ ) holds with e = Trgrr- ' 

Using the previous lemma we can state that an Hermitian operator is a per- 
turbation for a given seed S if and only if the following conditions are satisfied: 

(20) Supp(e) C Supp(S) 

(21) Tr[e7;^^]=0 Vmg5, Vi,j = l,...,m^ 



u 



(22) [e,Uh]=0 V/igG, 
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(conditions ( pi| ) and ( P2| ) follow directly from the normalization constraints (16) 
and ©). 

This set of conditions leads to an interesting property of extremal seeds: 

Theorem 1. R is an extremal point of C if and only if for any C G C one has 

(23) Supp(C) C Supp(S) =^ C-S. 

Proof. To prove necessity it is sufficient to define 6 = S — C and note that it is a 
perturbation of S. In fact, 8 is in the commutant of R(Go), Supp(8) C Supp(S), 
and Tr[8rj'^] =0 V/i G 5, Vi,j = l,...,r7i^. But, since S is extremal, then O 
must be zero. 

Viceversa , assume (P3). If 9 is a perturbation for S, then there exists some 



t ^ such that ^ = ^ + i8 e C. But a perturbation must satisfy (19), then 
Supp(C) C Supp(S). Using (JH) is then clear that 9 = ^-^(C - S) = 0. ■ 



The proposition tells us that extremal seeds have " minimal support" , in the 
sense that there is no element C G C with Supp(C) C Supp(S) which is different 
from S. 



Theorem 2. Let be ^ £ C. Write ^ in the form (\l'^). Then an operator Q is a 
perturbation of S if and only if 

(24) Tvi&T^^f] - V/x e 5, Vi, j = 1, . . . , m;, 
and 9 can be written as follows 

(25) 9 - (B.eso (^ ® XlA.X,) , 
with Xi, £ B{^^) and A,y E B{Rng{X^)) Hermitian V;/ G Sq. 

Proof. Suppose 9 is a perturbation. Condition ( pi|) is the same as (24). Due to 
condition (p2[), 9 must be an Hermitian operator in the commutant of R(Go), then 
we can write it in the block form 9 = (B^eSaili, (Xi Oi,), with each O^ G S(^^) 
Hermitian. Moreover, condition (|2y) along with ( [l7|) imply that each operator 
Oy must have Supp(0!y) C Supp{XlXi,) = Supp{Xi,). Using the singular value 

decomposition X^ = J2l=i ^^ l^i /(^i"! ({bDI s-'^d (l^i ) ^^^ orthonormal bases 
for Supp(Xi,) and Rng(X^) respectively) one can see that any Hermitian operator 
Oy with Supp(Oi,) C Supp(X^) admit the decomposition O^, = X^A^X,,, with A^, 



Hermitian operator in B(Rng{X^)). Conversely, if both conditions (^4|) and (25) 
hold, then conditions (EQ-E21) are obviously fulfilled. ■ 



Theorem 3. Let be P^, the projection operator onto the subspace J^ ^v Q J^ 
corresponding to the class v E Sq. An operator S € C written in the form S = 
©yg5p(/i, (g) X^X^) is extremal if and only if 

(26) e,e5„S(Rng(X,)) = Span{4^) \ ^^ E S, z, j = 1, . . . ,m4, 

where 

^ ■ ■ ^ Xl 



F-f = ©.e5o X, Tr,^^ 



P„T^'"^P,. 
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Proof. Using the characterization of Theorem g, we know that S is extremal if 
and only if for any operator 8 satisfying ( p4| ) and ( P5| ) one has = 0. Let's take 
Q in the form (Eq), and rewrite the direct sum as an ordinary sum 

(27) 6 = ^ P, (/, §5 XlA.X,) P., 

using the projectors P^ onto Jifl (g) M^. Using invariance of trace under cyclic 
permutations, we can write 



Tr 

(28) 



i^eSo 



e Tlf - 51 Tr {I. ^ A,){h ® X,)P,Tlfp,{h ® Xl) 



= ^ Tr \A, X,TT,^^[P,Tlfp,]Xl 



Define the space M = ®i,gSoRng(X^) and denote as (BveSo^^i^'^&iXv)) the hnear 
space of operators acting on ^ which are block diagonal on the subspaces Rng(X^), 
J^ G So- Then, the extremality condition for S becomes: for any Hermitian operator 
A e ®^eSol3iRng{X„)) one has 



(29) Tr 



AP/r^ 



V/iG S, Vi,j = l,...,m„ =^ ^ = 0. 



In terms of the Hilbert-Schmidt product {A, B) = Tr[Atp] this condition says that 
the unique Hermitian operator A G 0,ygSoS(Rng(Xy)) which is orthogonal to the 
whole set of operators F = {P,^^ | ^ G S, z, j = 1, . . . ,mfj} is the null operator. 
Orthogonality to the set F is equivalent to orthogonality to the set of Hermitian 
operators F' = {(P,f ' + PJf^) , z(P//' ~ P,'f') | A* G S, i,j = l,...m^}. Such 
orthogonality holds if and only if F' is a spanning set for the real space of Hermitian 
operators in 0,yg5f,S(Rng(X,y). Nevertheless, using the Cartesian decomposition we 
see that any complex block operator O G 0,yg5„S(Rng(X,y)) can be written as sum 
of two Hermitian ones, whence the extremality condition is equivalent to Span(F') — 
©j^gSo;B(Rng(X^). Finally, the observation Span(F') = Span(F) completes the proof. 
■ 

Notice that for trivial stability group Go — {e} (e denotes the identity element), 
we recover the characterization of H: there, one has indeed a single equivalence 
class V in So with one-dimensional representation space j?^, so that the whole 
Hilbert space M' is isomorphic to the multiplicity space ^p and the extremality 
condition dH) reduces to Span{Xr//'x^ | M G S, i, j = 1, . . . ,m^} = B{Rng{X)). 

Corollary 1. Any rank-one seed is extremal. 

Proof. Let be S a rank-one seed. In this case there is only one class i/q in the 
decomposition l ^u\ ) of S (otherwise S could not have unit rank), and the space 
B{Rng{Xi,g)) to be spanned is one dimensional, whence the condition ( p6| ) is always 
satisfied. ■ 

An alternative proof of Corollary ^ follows by observing that any rank-one ele- 
ment of the cone D of positive semidefinite operators is necessarily extremal for such 
cone: since the convex set C is a subset of D, a rank-one seed S G C is necessarily 
an extreme point of C. 
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Corollary 2. Let R Cz C be an extremal seed and write it in the form S = 

QueSoili' '^ XlX^). Define r^ = rank(X^). Then 

(30) T.-l^lK 

Proof. This relation follows directly from the extremality condition by noting that 
the left hand side is the dimension of the complex linear space of block operators 
®y(zSf^B[R'ng[X,y)), while the right hand side is the cardinality of the spanning set 

F = {^^^ ImGS, z,j-l,...,m4.B 

In Section g we will see an explicit example of extremal POVM which achieves 
this bound. 

5. Extremal POVM's and optimization problems 

A crucial step in a quantum estimation approach is the optimization of the es- 
timation strategy for a given figure of merit. This consists in finding the POVM 
which maximizes some linear (more generally concave) functional J- — e. g. the av- 
erage fidelity of the estimated state with the true one. Then, the convex structure 
of the set of POVM's plays a fundamental role in this problem, since, due to con- 
cavity of T ^ one can restrict the optimization procedure to the extremal POVM's 
only. 

In the covariant case, the problem resorts to optimize the state estimation in the 
orbit {UgpU^ I g E G} ~ G/Gq of a given state p under the action of a group G, Gq 
being the stability group of p. The optimization typically is the maximization of a 
linear functional corresponding to the average value of a positive function f{x, a;*), 
where the average is taken over all the couples (x, x^) of measured and true values 
x,x^ g X = G/Go, respectively. The joint probability density p{x, a;*) is connected 
to the conditional density p{x\x^,) given by the Born rule via Bayes, assuming an 
a priori probability distribution of the true value x^. In the covariant problem 
the function / enjoys the invariance property f{gx,gx^) = /(a;, a;*) Vg € G, and 
is taken as a decreasing function of the distance |a; — a;*| of the measured value x 
from the true one a:* . In the case of compact G one can assume a uniform a priori 
distribution for a:* values, so that the functional corresponding to the average can 
be written as follows 

(31) Tp[E] = f dg f dg, f{gxo,g.xo) TiiUg,pUlUgEUl] 

Jg Jg 

(32) = f dg fixo.gxo) Tr[UgpUlE], 

Jg 

where xq is the equivalence class containing the identity. In the following, we will 
consider as the prototype optimization problem the maximization of the likelihood 
functional |, | 

(33) Cp[E] = Tr[pS], 

corresponding to the choice /(a;, a;*) = 6{x — a;*) in Eq.(|3|). Maximizing Cp[E] 
means maximizing the probability density that the measured value x coincides 
with the true value a;*. For such estimation strategy the optimization problem has 
a remarkably simple form, enabling a general treatment for a large class of group 
representations |13]. Moreover, the solution of the maximum likelihood is formally 
equivalent to the solution of any optimization problem with a positive (which, a 
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part from an additive constant, means bounded from below) summable function 
f{x,x<,). Indeed, we can define the map 

(34) M{p) = fc-i / dg f{xo,gxo) UgpU^ 

where k = Jg'^9 f{xo,gxo). This map is completely positive, unital and trace 
preserving, and, in particular, M[p\ is a state. With this definition, we have 

(35) J^p[S] = fc £^(p)[S], 

whence the maximization of Tp is equivalent to the maximization of the likelihood 
for the transformed state M{p). 

Essentially all optimal covariant measurements known in the literature are rep- 
resented by rank-one operators. The rank-one assumption often provides a useful 
instrument for simplifying calculations. Nevertheless, as we will show in the follow- 
ing, the occurrence of POVM's with rank grater than one is unavoidable in some 
relevant situations. 

Proposition 1. For any S e C, 

(36) rankfSl > max ( —r- 

Proof. Let's decompose J^ into irreducible subspaces for the representation R(G) 
of G as follows 

(37) ■y^ = ®^.^s®Z\'^}''^■ 

Take an orthonormal basis B^-'^-^ — {|(/z,i),n) | n = l,...,(ip} for each subspace 
J^}^^ in such a way that \{iJi,i),n) = T^f\{ii,j),n) for any n, T^^f : Jf^ ^ Jf, 
being the invariant isomorphism which intertwines the equivalent representations 
{p,i) and {fJ-,j). Diagonalize S as 

rank(H) 

(38) S= E l^^)('?'=l 

fe=i 

and write 

(39) W-EEE^(M.).nl(A^'*)''^)- 

fi^S i—l n—1 

Since iml^ijlVk) = J2nLicf^A),n^tfi.j),n^ ^^^ normalization constraints ((l6|) be- 
come 

rank(H) d^ 

(40) E "^ 4^.i),n4^^J),n=d^.S,J■ 
k=l ri=l 

This relation implies that for any /i G S the vectors {c(^ni\ \ i ~ 1, . . . , ni^} defined 

by {cr^,i))k,n = cf j\ „ are orthogonal: since they are ?n^ orthogonal vectors in a 

linear space whose dimension is d^ x rank(S), it follows that tti^ < d^ x rank(S), 

hence rank(S) > ^ V/^ G S. ■ 

Summarizing, every times m^^ > d^ for some class /x G S, a covariant POVM 
cannot be represented by a rank-one seed. 
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The previous proposition exhibits a structural reason for which, in the presence 
of equivalent representations, the set C of covariant seeds may contain only elements 
with rank greater then one. On the other hand, in the following we will discuss the 
occurrence of covariant POVM's with rank greater than one in explicit optimization 
problems, independently of the presence of equivalent representations. 

Proposition 2. Let he S an extremal point of C. Denote by P the projector onto 
Supp(S), and let r = rank(P). Then S is the unique seed which maximizes the 
likelihood for the state p — —. 

Proof. 

First, we need to prove that S commutes with the representation R(Ho) = 
{Uk I fc € Hq}, where Hq is the stability group of p, defined by [p, Uk] =0 Vfc G 
Hq. Define the group average 

Since R(Ho) is the stability group of the projector onto Supp(S), clearly Supp(S) 
is invariant under R(Ho), whence £, satisfies Supp(^) C Supp(S). Moreover, using 
the invariance of the Haar measure it is easy to see that ^ commutes with R(Ho). 
Finally, ^ is an element of C. In fact, it is positive semidefinite, satisfies (|l6|) and 
commutes with R(Go) — the stability group of S — which is by definition a subset of 
R(Ho). Since S is extremal, using Theorem [l| we can conclude that S = ^, whence 
S commutes with R(Ho). 

Let's prove now optimality. For any arbitrary seed C G C, the following bound 
holds: 

,42, ,,M.T.Ki^a£fl<^^^^=«. 

r r r 

where the last equality follows from the normalization constraints dlq). Clearly S 
achieves the bound, whence it is optimal. Notice that the inequality Tr[PC] < Tr[C] 
becomes equality if and only of Supp(C) C Supp(S), then using Theorem ll| we can 
see that S represents the unique optimal POVM. ■ 

Consider now a density matrix a with support in the orthogonal complement of 
Supp(S), and consider the randomization 

P 

(43) p={l-a) \-aa, 

r 

with < a < 1. In the following we prove that, for sufficiently small a > 0, S 

is still optimal for the maximum likelihood strategy. In other words, the extremal 

POVM represented by S is stable under randomization, and the same measuring 

apparatus can be used for a larger class of mixed states. 

Proposition 3. Consider the randomized state p in (um and denote by q the 
maximum eigenvalue of a. If a < i^r- r then S is the unique seed which maximizes 
the likelihood for the state p. 



Proof. First, notice that ^ commutes with the representation R(Ho) of the stability 
group of p. This follows from the observation that the condition a < .]; _ implies 
that ^—^ is strictly the largest eigenvalue of p. Then, P is the projector on the 
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eigenspace with maximum eigenvalue of p, while, for any h E G, Ph = UhPUj^ is the 
projector on the eigenspace with maximum eigenvalue of ph = UhpU^- If /i G Hq 
then it must be pu = p, and, necessarily, Ph = P. Therefore Hq is a subgroup of 
the stability group of P. But S commutes with the representation of the stability 
group of P, as proven in Proposition 0, then it commutes also with R(Ho). 

Now we prove optimality of S. Let's denote by Q the projection onto Supp((T). 
The following bound holds for any C g C: 

(44) Cp[C] = ii^Tr[PC]+«TrK] 

(45) < ^"~"^ Tr[PC]+agTr[QC] 





r 




(1 


- 


a) 




r 




(1 


- 


a) 




r 




(1 


- 


a) 



(46) < ^— ^ Tr[(P + Q)C] 

(47) < ^^ — ^ Tr[C] ^ ^^ ~ "^ dim(,jf ). 

r r 

This bound is achieved by S, proving its optimality. Notice that S is the unique 
optimal seed. In fact, equality in (|46|) is attained if and only if Tr[(5C] = 0, namely 
when Supp((3) C Ker(C), while in ( ]47|) equality is attained if and only if Supp(C) C 
Supp(P) ® Supp(Q). Therefore the bound is achieved if and only if Supp(C) ^ 
Supp(P) = Supp(S), implying C = S. ■ 

6. Examples 
6.1. Extremal POVM's with a non trivial stability group. 

6.1.1. Consider the group of rotations, represented in a (2j + l)-dimensional 
Hilbert space J^j by the irreducible representation Rn.tp = e"^"J, where ip is an 
angle, n is a unit- vector, and j = {jx,jy,jz) is the angular momentum operator. 
In this case a covariant estimation in the orbit of a pure state \^) generally may 
involve a nontrivial stability group. This is actually the case when |^) = \jm)no, 
is an eigenvector of no ■ j for some unit vector Uq. Clearly in such case the sta- 
bility group Go consists of rotations around no, and the state estimation in the 
orbit reduces to the estimation of a rotated direction n'. The same situation arises 
for any state p mixture of eigenvectors of no ■ j. Without loss of generality, let's 
take no as the direction of the z-axis, and write p = X]m=-7"P"il.?'^)0'^l with 
Pm ^ Wm. Let's denote by P the projector onto Supp(/9), and take m such that 
Pfn = max„{p„}. Then, since 

TrK] < Prn Tr[PS] < p^ Tr[S] = p™(2j + 1), 

one has that S = (2j + l)\jrh){jrh\ is the optimal POVM. Notice that such POVM 
commutes with the stability group R(Go) and is extremal, as a consequence of 
Corollary ^. 

6.1.2. Consider the group §U((i) of unitary d x d matrices with unit determinant, 
acting on the space Jf^ ^ C*. It is easy to see that each vector \^) S J^ has a 
nontrivial stability group Go = U(d — 1). In fact, by introducing an orthonormal 
basis B_L = {|n) |n = l,...,(i— 1} for the orthogonal complement J^^ of the line 
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Spanjl?/;)}, and the basis B = |-(/;) U B_l for ^, the stabihty group Go consists on 
matrices of the form 

(48) Uu-- ^'^'^ 







Vh 



where ojh G C, \uJh\ = 1, and Vh is a unitary [d— 1) x (d— 1) matrix with Det(Vh) = 
w^. Let's consider now the tensor representation R(G) ~ {Uf^ \ Ug G §U((i)} on the 
space M"^"^ . This representation has two irreducible subspaces, the symmetric and 
the antisymmetric ones M+ and J^- , with dimensions d+ = ^ ^ ' and d_ = ~ ' 
respectively. Denote by P+ and P_ the projectors on J^ and J^. Let's apply the 



representation R(G) on the state |V')® G -J^^ . Clearly the stability group is the 
same Gq as before, and it is represented by R(Go) = {Uf^ \ h G Go}. It is easy to 
see that R(Go) contains five irreducible components, carried by the subspaces Mi = 
Span{|i/>)®2}, ^2 = Span{|i/;)}0^^ , =^3 = ^^(g)Span{|V)}, ^4 = P+[M'^ ®2)^ 
and ^ = P_(^^ ®^). Notice that ^ and J^ carry equivalent representations, 
corresponding to a two dimensional multiplicity space. An example of extremal 
POVM is given by 

s-^l^)(^re^P-OP-, 

where Q is the projection on ,^^ ^"^ . Since the two summands are proportional to 
|^)(i/'|®^ and P-QP-, which are the projectors on Jtfi and J^5 respectively, then S 
belongs to the commutant of R(Go) = {U^'^ \ h £ Go}. Notice that the subspaces 
■J'^i and ._^ have multiplicities mi ~ 7715 — 1, corresponding to one-dimensional 
multiplicity spaces ^i = ^5 = C (whence the partial traces over J^i^^ will be c- 
numbers). Moreover, using the fact that Tr,^JP+] = 1, Tr^jP_] = 0, Trj)^g[P+] = 
0, Tr^5[P_] = 2 ^^^ '^^^ check extremality using the condition (|2q). Let's 

observe that in this example we have ri — r^ ~ \ and 777+ = 777_ = 1, where ri and 
rs are defined as in Corollary 0, while 777-|_ and 777_ are the multiplicities of the two 
irreducible representations of R(G). Then the bound of ( pO| ) is saturated. Finally, 
we remark that this POVM is optimal for discriminating states in the orbit of |-0)®^ 
(13), in the orbit of p = i (|i/')(V'|®^ + P-QP-) where r = 1 + ('^-^K'^-^) ^^^^^^gg ^f 
Proposition y, and also in the orbit of any randomization p' ~ {1~ a)p + aa where 
(T is density matrix with Supp((T) C Ker(P), and a < j-^, because of Proposition 

I 

6.2. Extremal POVM's with rank greater than one. 

6.2.1. Consider the Abelian group G = U(l) of phase shifts, acting in the space 
^ = C* by the representation R(G) = {U{(p) = exp{iipN} \ Lp G [— 7r,7r]}, 
where the generator N is given by A^ = Sn=o "■ l"')(^l ^'-'^ some orthonormal basis 
{\n) I 77 = 0, 1, . . . , d — 1}. The stability group Gq may be either the whole U(l) 
(for p diagonal on the eigenstates of the generator), or a discrete subgroup Go = Zj. 
for some integer fc, including the case fc = 1 of trivial stability group. We exclude 
the degenerate case Go = U(l) of shift invariant states. The parameter space 
X = U(l)/Zfc will be a circle, parametrized by an angle 9 G [— 7r,7r], and the action 
of a group element g{(p) G G on an element G X will be given by g{i^) 9^9 + kip. 
Due to constraint (|6|), a seed S is represented in the eigenbasis of the gener- 
ator by a correlation matrix, namely by a positive semidefinite matrix with unit 
diagonal entries. Viceversa, any correlation matrix corresponds to a seed in the 
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case of trivial stability group Gq . In [Q one can find a constructive method which 
provides extremal correlation matrices with rank r > 1: here we show that any 
of such matrices can be viewed as the optimal seed for the estimation problem 
in the orbit of a particular state. Let us choose as optimality criterion the max- 
imization of the average value of a positive summable function / : X x X ^ M+ 
depending only on the difference 6 — 9^ between the measured and the true value. 
Suppose p a state with stability group Go = l^k- As we noted at the begin- 
ning of section pi the maximization of Fp[S] — the average value of f{9 — 6**) — 
corresponds to the maximization of the likelihood /!^(p)[S] for the transformed 
state Mip) = /o"^ /^^ ^f{-kLp)U^pUl (from Eq. (|3|)). Notice that the map M 
is trivially covariant — i.e. M.{U^pU^) = U(j,M{p)U^ — since the group is abelian. 
For simplicity here we require that the map A4 is invertible, whence also Al^^ is co- 
variant and trace- preserving (but generally not positive). Covariance of A^ implies 
that the stability group of M {p) contains the stability group of p, and covariance 
of M^^ implies the reverse inclusion, whence the stability group is not changed by 
the maps. 

Let's take now an extremal correlation matrix S with rank(S) ~ r > 1 and 
denote by P the projector onto Rng(S). Using Proposition pi we can see that S 
commutes with the representation R(Ho), where Hq is the stability group of P. 
Call A the modulus of the minimum eigenvalue of A^^^(— ), then 

is a density operator. Notice that the stability group Gq of p is the same stability 
group of A4^^{P), which coincides with Hq, the stability group of P. Therefore S 
commutes with the representation R(Go). It is easy to show that S is the unique 
seed commuting with R(Go) which is also optimal for the estimation of states in 
the orbit of p. In fact, for any C in the convex set C of the seeds with stability 
group Go, we have 

F,[C] = /oTr[CM(p)] = /of^^Tr[C] + ^7^^Tr[CP] 



/ , :d\ fl + rX 



1 + d\ '^' r(l + dX) 

1 + dx) 

This bound is achieved choosing C = S, moreover, as in Proposition y, we can 
observe that the functional Tr[CP] with ^ G C is maximum if and only if C = S, 
then the maximum is unique. 

6.2.2. We provide now an example with a non-compact group represented in an 
infinite dimensional Hilbert space. This example is out of the general treatment of 
the present paper — which considers only finite dimensions — and is given only with 
the purpose of showing that our results could be generalized to infinite dimensions, 
however at the price of much more technical proofs. 

Take J^ as the Fock space, and consider the projective representationon Jif of 
the group of translations on the complex plane C in terms of the Weyl-Heisenberg 
operators R(G) = {D{<^) = e"" """ | a £ C}, where [a,a'^] — 1. Here we will con- 
sider the 2-fold tensor representation {£)(a)®^ | a g C} on J^^'^. Using the unitary 
operator V = gf(a.i4-o.\a2) ^ ^^^ ^^^ ^j-^^g D{a)^^ = V{D{V2a) (g) I)V'' and see 
that the irreducible subspaces of this representation are J^n = V{^ ® Span(|(/)„)), 
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{{(f'n) I n — 1,2, . . .00} any orthonormal basis for J^. All these subspaces carry 
equivalent representations, the isomorphism between J^m and J^n being 

(49) r„„ = i/(/®|(/.„)(0„|)yt. 



In terms of these isomorphisms, the normalization constraints (Iq) for a seed oper- 
ator become Q 

(50) Tr[r,„„C] - 2(5„„, 

Notice that the number 2 in this formula has nothing to do with the dimension 
of J^n which is infinite: in the non-compact case the dimensions are replaced by 
positive numbers depending only on the equivalence class of representations. In 
principle, since the space J^®^ is infinite dimensional, there is the possibility of 
extremal covariant POVM's with an infinite rank. Actually we can provide the 
remarkable example 

(51) E = 2Vi\0){0\(g>I)V\ 

where |0) is the vacuum state of the Fock basis {\m) \ a^aJTO) ~ m\m)}. The 
corresponding POVM can be realized by averaging the outcomes of two independent 
measurement with Si = |0)(0|(Xi/ and S2 — /(81 10)(0| |13], which in quantum optics 
correspond to two heterodyne measurements |l4] . 

We can observe that S is maximizes the likelihood functional for any state of 
the form p = F(|0)(0| ® a-)V\ where a ~ X)^oP"l<^")('^nl' ^^ ^ mixed state with 
Pn > Vn. In fact, for any seed C, one has the bound 

(52) 

OC 

TT[V{\0){0\^a)V^C] =Y,PnTr[V (|0)(0| ® |(^„)(0n|) V^ Q 

00 00 

<^p„Tr[T/(/® |0„)(<^„|)0C] = 5IP"Tr[r„„C] = 2, 

n— n 

and since S achieves the bound (p2|), it is optimal. Moreover S is the unique 
optimal seed. In fact, the equality in ( |52| ) is achieved if and only if Tr[F(|0)(0| (g) 
|0n)(0n|)l^^C] = Tr[V{I IS) |0n)('/'n|)V"^C] for any n: by expanding the identity on 
the Fock basis, the positivity of C implies {m\{(f>n\V''CV\m)\(f)n) — for any m 7^ 
0. Hence the unique nonzero diagonal elements of ( are on the vectors V^|0)|(/>„). 
On the other hand, the positivity of C along with the normalization constraint 
Tr[rm„C] — Vm 7^ n imply that all the off diagonal elements of C are zero. 
Hence C = 21^ ^^^^(|O)(O|®|(/)„)(0„|)Ft = 2F(|0)(0| (^/)T/t = S. The fact that S 
is the unique optimal seed ensures that it is also extremal, otherwise there would 
be two different seeds which are equally optimal. Notice that S is extremal also 
according to our characterization (|2q ) . 
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